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Abstract 

We analyze the stress-energy tensor, and the resulting energy conditions, for a scalar field with 
general curvature coupling, outside a perfectly reflecting sphere with Dirichlet boundary conditions. 
For conformal coupling we find that the null energy condition is always obeyed, and therefore the 
averaged null energy condition (ANEC) is also obeyed. Since the ANEC is independent of curvature 
coupling, we conclude that the ANEC is obeyed for scalar fields with any curvature coupling in 
this situation. We also show explicitly how the spherical case goes over to that of a flat plate as 
C/3 \ one approaches the sphere. 
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I. INTRODUCTION 



From our everyday experience we expect the energy density of any substance to be a 
positive quantity. However, quantum field theory allows the existence of states which have 
negative energy densities, as first shown theoretically by Hendrik Casimir in 1948 1]. Casimir 
found that the energy density between two neutral conducting plates was negative and 
resulted in an attractive force between the plates. This attractive "Casimir" force has only 
recently been accurately measured 0. 

The reality of negative energy densities has opened the possibility of exotic physical 
situations in general relativity, including wormholes and time travel. With no restriction on 
the stress-energy tensor, one can simply write down any desired space-time geometry and 
solve Einstein's equations in reverse to find the matter stress-energy that will produce that 
geometry. To prevent exotic phenomena requires energy conditions that all matter must 
obey. 

The simplest of these is the weak energy condition (WEC), which states that no observer 
sees negative energy density, 

T^WVtO, (1) 

where T^ v is the stress-energy tensor and is any timelike vector. If instead we consider 
null in Eq. ([Tfl. we get the null energy condition (NEC). If the WEC is obeyed at some 
location for all V, then by continuity the NEC must also be obeyed. These conditions, 
however, are violated by the Casimir effect between parallel plates. 

The averaged weak energy condition (AWEC) allows the WEC to be violated locally, but 
holds that when one integrates the WEC over a complete geodesic with tangent vector V^, 
the result must be non-negative, 

roc 

dxT^VV >0. (2) 

When is null, Eq. (J2J) represents the averaged null energy condition (ANEC). The AWEC 
is also violated by the Casimir effect, but the ANEC is obeyed, because the geodesic cannot 
pass through the plates (considered as ideal boundaries) and so must run parallel to them, 
in which case all contributions cancel. 

The ANEC is the most important condition for our purposes, since its validity is sufficient 
to rule out many exotic phenomena, such as traversable wormholes j^, super luminal travel 
0], and closed timelike curves [15], and to prove singularity theorems SS] 1 - No ANEC 
violations are known for a scalar field in flat space quantum field theory. In [9J], we computed 
the stress-energy tensor outside a Dirichlet sphere and found that ANEC was always obeyed. 
Here, we extend this work to arbitrary curvature coupling. 

The stress energy tensor for an arbitrarily coupled scalar field outside a sphere was found 
by Saharan [10]. We reproduce these results using a similar technique, based on the calcu- 



lational framework developed in ll(. In Sees. HP and UTT1 we obtain a general expression for 
the vacuum stress-energy tensor outside a spherically symmetric scattering center. In Sec. 
UVI we specialize to perfectly reflecting boundary conditions on a sphere. In Sec. |VJ we show 
how our results give analytic expressions for the stress-energy tensor in the large r limit. 



1 For superluminal travel, we integrate only over the path to be traveled, and for singularity theorems only 
over the part of the geodesic to the future of a trapped surface. 
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In Sec. IVII we calculate the NEC numerically and find that it is always obeyed for confor- 
mal coupling. Since the NEC is obeyed for conformal coupling at every point outside the 
sphere, we know that the ANEC is obeyed too. Furthermore, the ANEC is coupling constant 
independent, so we can conclude that the ANEC is obeyed for all scalar field couplings. 

In Sec. IV 111 we show how our results approach the usual stress-energy tensor for a flat 
Dirichlet plate for points very close to the sphere. We follow with conclusions in Sec. IV 1111 



II. STRESS-ENERGY TENSOR FOR GENERAL COUPLING 



The stress-energy tensor for a massless real scalar field is 



- d,dM 2 )] , (3) 

< "^Planck- For 



where £ is the curvature coupling, and we are working in the limit where 
conformal coupling in 3+1 dimensions, £ = 1/6. 

If we consider a static, spherically symmetric system, the stress tensor can be written in 
terms of the energy density p, the radial pressure p r and the tangential pressure p±. For 



example, on the x axis we have T 00 
vanishing. 

From Eq. (j3J), these components are 



Pi T xx p r , X 1 , 



yy 



p±, with other components 





1 " 


p = 


2 . 




1 " 


Pr = 


2 . 




1 " 


P± = 


2 . 



2 + (d r 4) 2 + 2(d ± <p) 2 
2 + (d r 4) 2 - 2(«9 ± 0) 2 j +£[2dl( 
2 -(^0) 2 l+e[9i(0 2 )+9 2 (0 2 



+ 2d 2 



(4) 
(5) 
(6) 



where d± denotes the derivative along a tangential line. (By rotational symmetry the choice 
of line will not matter.) The vacuum expectation value (<j) 2 ) depends only on r, since 
the system is static and rotationally symmetric. Thus (d 2 ((p 2 )) = 9 2 (0 2 ) = 0. The second 
derivative in the tangential direction does not vanish, however. There is still a term dj_ ((f) 2 ) = 
(l/r)<9 r (0 2 ). Thus 

' U 2 + («9 r 0) 2 + 2(«9 ± 0) 2 



P 

Pr 
P± 



2 + (d r 4>) 2 - 2(d ± C 



2^ 



r 



2 - (drCj))' 



-d r 
r 



dm 

d r 

dlun 



-Or 



(7) 
(8) 
(9) 



If we are interested only in the null energy condition, then the terms proportional to the 
metric do not contribute, and we have 



T„„V»V V = (V a d a 4>) 2 - ZVVd&tt 2 ) . 



(10) 



If V 



'l,v) with Ivl = 1, then 



tf + vtidr^ + vlid^) 2 -^ 



V 2 r d 2 M 2 ) + -9r 



where v± is the magnitude of the velocity in the tangential direction, v\ + v. 



(11) 
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III. VACUUM EXPECTATION VALUES 



In order to calculate the vacuum expectation value of the stress-energy tensor, it is 
sufficient to compute the expectation values of (ft 2 , 2 , (<9 r </>) 2 , and (d±4>) 2 . 
We can write the quantum field <fi in spherical modes 

#r,n,f) = X;£/ ^^=^(r)^(^)*e^aL m + c.c. (12) 

— — /n \ 1X0J 

where ui = k since the field is massless. 

The wavefunctions ipi{r) satisfy the time-independent radial Schrodinger equation 



d 2 2d 



dr 2 r dr r 2 



+ -^)rt(r) = k 2 rt(r). (13) 



normalized 

poo 

k 2 r 2 ^ k {ry^ k ,{r)dr='-5(k-k') (14) 



00 — 

7T 



2 

so that the free wave function is just 4>k( r ) = tt{k r ) where j is the spherical Bessel function 
j t (z) = y/ir/(2z)Jt +1 / 2 (z). 

From Eq. (|12|) we can compute the vacuum expectation value of <ft 2 . There will be no 
dependence on the angular position, so we can sum over the spherical harmonics using 



m=-t 

to get 



^ 2 ) = E^ r dk^- (\4(r)\ 2 - h(krf) (16) 

where we have subtracted the free wavefunctions. No other counterterms are necessary, 
because we are considering locations in empty space without any potential. 

The vacuum expectation values for the time and radial derivatives are straightforward, 

(<P 2 ) = J2 2 ~T^ rdkk>u,(M(r)\>-Mkr)*) , (17) 
t ^° 

_ Of 1 1 foo 7,2 

<^) 2 > = £ 4^T / dk ~ (I^WI 2 - (drUkr)) 2 ) . (18) 

To compute the azimuthal term we need the derivative of the spherical harmonics. With- 
out loss of generality we can work in the equatorial plane where we can take d±Yg m (jr/2, 0) = 
(l/r)d<f,Ye m (7r/2,<j)) = (im/r)Y em (7c/2,(j)). Using the formula 

±^Y lM )^ eil+ T l+1) (19) 
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we get 

roo 1 2 



((MVE ^W^ r^^(l^(r)| 2 -^(Ar) a ) . (20) 
Outside any spherically symmetric potential, the wave functions are given by 

^(r) = X - Y ih hf (At) + hf (At)] (21) 



where is the scattering phase shift in the quantum mechanical problem with the same 

'2) 

1/2 



potential, and the spherical Hankel functions are hf (z) = y/nz] (2z)H^\ 2 (z). Thus 



l^f(r)i 2 -h{krf = \ [(e 2 ^ - l) h^\krf + (e~^ - l) h^\krf] . (22) 



From the appendix of [9J] it is easy to show that when we use Eq. (|22jl in Eqs. 
we can drop the second term of Eq. (|22j) and extend the range of integration over k to — oo, 
with the understanding that k is to be taken above any branch cut on the negative real axis. 

Then, following the methods used in and [§] (a similar procedure was used previ- 
ously in [13]), we convert each expression to a contour integral which we close in the up- 
per half plane. The only contribution to the integral comes from the branch cut along 
the imaginary k axis. To the right oj = \fk^ = k, but to the left oo = —k. We 
then use hf\ix) = —(2/7r)i~ e ke(x), where kg is the modified spherical Bessel function, 
ki(z) = a/tt/ {2z)Kf + i/ 2 {z). With k = in, we obtain 

{<P 2 ) = - r dK K - !) (23) 



(0 2 ) = E(-)'^ 1 r *« - 3 - 1) **w a ( 24 ) 

<(^) 2 > = - jf rfK - W 2 ( 25 ) 

((M) 8 ) = " EH ^ + 4 ^f +1) r d« - - 1) W . (26) 

IV. PERFECTLY REFLECTING BOUNDARY CONDITIONS 

For a hard sphere of radius a, with Dirichlet boundary conditions, 



e 



2is e = hf\ka) 
hf\ka) 



(27) 



so that 6(a) = 0. Thus 



and 



e 2i5e - 1 = - 7" v ^ y (28) 



_ ! = (_)^^W (29) 
ke{Ka) 



since je(ix) = i e ie(x) where ig(z) = a/ 71 "/ (2z)Ie + i/ 2 \ 
Thus 



z ■ 



= _ <( < + i )(2< +1) r^M^ 

^T -1 47TT 2 y 



(30) 
(31) 
(32) 
(33) 



in agreement with (Io| . 



agieeuieiiL wim [iuj. 

From these equations we can compute any combination of stress-energy tensor compo- 
nents. For V null, from Eq. fTTf . 

,2 



z z 

V -±-cP r {h{Kr) 2 )+^d r {h{Kr 
k z K z r 



V. LARGE DISTANCE BEHAVIOR 

Far away from the sphere, in the large r regime, the stress-energy tensor is dominated by 
the I = partial wave contribution. For this case the integrals in Eqs. (I30H33|) can be done 
analytically. We find 

„2 



l.2\ 



((dA) 2 ) 
{{d^f) = 



ln(l — a/r) 
g 7r 2 r 2 

2r — a 

32vr 2 (r - a) 2 r 4 " 16vr 2 r 5 ' 32vr 2 r 6 
5 — 6r + 4(r — a) 2 ln(l — a/r) 5a 



87r 2 r 3 167r 2 r 4 

a 3a 2 
+ 



9a 2 



32vr 2 (r - a) 2 r 4 



167T 2 r 5 327r 2 r 6 



(35) 
(36) 

(37) 
(38) 



where the approximate forms are those necessary to compute components to 0(r~ 6 ). 

Substituting these large r approximations into Eqs. (J7HHJ), we find the following stress- 
energy components 



P 

Pr 
Pi 



3a 2 



37r 2 r 5 



32ttV 
3a 2 



8vr 2 r 5 327rV 
3a 3a 2 



+ 



167r 2 r 5 167r 2 r 6 



3a 3a 2 
+ 



47r 2 r 5 47r 2 r 6 



3 a 



.2 1 



+ 



9a 



+ 



(39) 
(40) 
(41) 
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The leading terms agree with [Io| . with the correction that the fraction in Eq. (4.17) of 
Tol | should be {D — 1)/D rather than (D — 2)/(D — 1), where D is the number of spatial 
dimensions [l2j. 

For conformal coupling all terms of order (1/r 5 ) cancel to give 



Using Eqs. 



cr 



P 



32vr 2 r 6 

,2 



cr 



Pr = 
P± 



96n 2 r ( 

,2 



a 



48vr 2 r 6 ' 

in Eq. |TT| we find a large r expression for the NEC, 
3a 2 f 5a 9a 2 



a 



167r 2 r 5 

,2 



+ 



v.,. 



327r 2 r 6 
15a 



+ 



- v r 

3a 2 
2vr 2 r 6 



+ 



167r 2 r 5 327r 2 r 
3a a 



4 7r 2 r 6 



For conformal coupling, 



327r 2 r 6 

so the NEC is obeyed far from the sphere. 



(42) 
(43) 
(44) 



(45) 



(46) 



VI. NUMERICAL RESULTS FOR THE NEC WITH CONFORMAL COUPLING 

We have computed the NEC contribution for the case of conformal coupling, £ = 1/6, 
by numerically performing the integral in equation Eq. (|34p. Fig. ^ shows the contributions 
from several partial waves for the case of radial motion. The case of tangential motion is 
very similar. The i = contribution is always positive, as shown in Eq. (|46|) . but for other 
£, the sign depends on position. Thus, unlike the electromagnetic field case [l3j], the NEC 
does not hold in each channel individually. 

Nevertheless, when we sum up all contributions, we find that the NEC is obeyed at any 
distance from the sphere, as shown in Fig. EJ Tangential motion yields a larger value for the 
NEC at any given point than radial, but both are always positive. Since the NEC is obeyed 
at each point, the ANEC is also obeyed. 

The term depending on £ in the NEC, Eq. (11 0|) . is a total derivative in the direction V, so 
it gives no contribution to the ANEC integral in Eq. (j2J). Thus the ANEC does not depend 
on curvature coupling, and since it holds for conformal coupling it will hold for any coupling, 
as we found numerically in the minimally-coupled case 0. 



VII. PERFECT MIRROR RESULT FROM SPHERICAL CASE 

In this section we will consider the leading-order behavior of the stress-energy tensor 
when we get very close to the sphere, the opposite limit to that of Sec. El In this case, there 
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FIG. 1: We plot here the null energy condition, at a distance r from the center of a sphere for 
motion in the radial direction for the i = (solid line), 1 = 2 (dashed line), I = 4 (dashed-dotted 
line), and t = 6 (dotted line) partial waves. 



will be a divergence due to summing over large numbers of partial wave modes, so we can 
use the large-order approximation for Bessel functions, as described in (lp| . 
The expectation value of 2 is given by Eq. ([30)1 . which we can write 

where u = £ + 1/2. 

Now we will consider the case where e = r/a-l<l. We let x = na/u and y = (l+e)x = 
nr/v to get 

<« = -E^f^^W («) 

Since we are summing over a large number of partial waves, we can use the large-order 
approximation to the Bessel functions. For z/> 1 Q, 

1 e vri 

h(vx) " ^(TW (49) 



7T 



K ^ X) *> ^(TW 1 (50) 
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1.2 1.4 1.6 1.8 



FIG. 2: The null energy condition at a distance r from the center of a sphere for motion in the 
azimuthal (solid line) and radial (dashed line) directions. 



where 



so 



7}{x) = y/i + x 2 + In 



x 



i + VTTx~ 2 



(0 2 



An 2 a 



-T 

■2 n 2 Z—/ 



V 



dx 



-2v(r)(y)-r)(x)) 



o vTT 



(51) 



(52) 



We have ignored higher order terms in e by replacing r with a in the prefactor and y with 
x in the square root. 

We expand to first order, 



dr) 



v(y) ~ v( x ) = tx~r = tVl + x 2 
dx 



so we have 



Thus 



dx 



-2ueVl+x 2 



o VTT 



K {2ve) 



^)--Ew^(2, e ). 

We approximate the sum by an integral, 

f°° v 1 



2 e 2 167T 2 (r — a) 



(53) 

(54) 
(55) 

(56) 
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From Eq. (|31|). the time derivative term differs from (0 2 ) by a factor — k 2 in the integrand. 
By the analysis above, 



z/ 3 r , x 2 



Using 



we find 



3 ) » V — / e ~Mn(y)-n(»)) (57) 



°°^^^e--^=^M (58) 







Vl + x 2 2z/e 



z/ 2 



(6 2 ) » y — ;— - KA2ve) « / aV = . K x (2ve\ = ^-r = ^ rj (59) 

Vr/ ^87r 2 a 4 e V ; 7 8vr 2 a 4 e V ; 32vr 2 a 4 e 4 32vr 2 (r - a) 4 V ; 

Since ~ l/z for small 2, the contribution from a single partial wave mode with 

£ C 1/e is proportional to £/e 2 . The sum over these modes thus diverges as £ 2 . It is 
truncated when £ ~ 1/e, so the final result is proportional to 1/e 4 . 

From Eq. (J33J), the tangential term differs from (0 2 ) by a factor £(£+ l)/(2r 2 ) « z/ 2 / (2a 2 ) 
since we are keeping only the highest power of v and lowest order in e. Thus we find 

<(^) 2 > « - / A/jTr^o^i/e) = -^m = ^ (60) 

J m z ar 327ra 4 e 4 327r^(r — a) 4 

The angular derivative brings in a factor of £, so this term has two more powers of I and 
correspondingly two fewer powers of e in the denominator. For I -C 1/e it goes as £ 3 In z/e, 
which is much smaller than the time-derivative term, even though the sum over £ is the 
same except for sign. 

From Eq. (|32|). the radial term differs from ((J) 2 ) by a factor —k'^nr) 2 / ki{nr) in the 
integrand. The derivative of the spherical Bessel function is 

£ £ + 1 

= ~2£~Tl ke+1 ^ ~ 2£~Tl ke ~ 1 ^ ^ 

while the derivative of the regular Bessel function is 

K' l (z) = - l -K t+l {z)- l -K t _ 1 {z) (62) 



but the difference in the coefficients does not matter for £ large. Thus 

k e (z) ~ K t (z) 

The derivative for large order is 

~(^ _l „2U/4 . /i _l ,,2 



(63) 



*IM » = (64) 

2u y y 




Thus we find 



k'JKr) 2 1 



z/ 



2 



1 + ^ = 1 + 7^ ( 65 ) 



kiinr) 2 y 2 K 2 r 2 

10 



Consequently 

({d r <py) « 2((9 ± 0) 2 ) - (<p 2 ) 

For a conformal field we will also need 



327r 2 (V - aY 



7r 2 (r — aY 



(66) 



(67) 



from Eq. ([Hfijl. but we will not need d r ((f> 2 ) because that does not diverge as (r — a) 4 . 
Putting everything together, the stress tensor components to order l/(r — a) 4 are 



P 

Pr 
P± 



<p 2 + (d r <p) 2 + 2(d ± <p) 2 

<p 2 + (d r 4>) 2 - 2(d ±( p) 2 
4?-{d r M +Zd 2 (<p 2 ) 



-eW) 

« 



l-6£ 



16vr 2 (r - aY 



6£ 



167r 2 fr — a) 



(68) 

(69) 
(70) 



For minimal coupling we get the standard result for a single, flat Dirichlet plate, and for 
conformal coupling, £ = 1/6, we also get the standard vanishing result for the flat plate. 
Therefore, by linearity, the result is correct for general coupling. 



VIII. CONCLUSIONS 

We have studied the problem of a real, massless scalar field outside a reflecting sphere, 
in the calculational framework developed in We obtained expressions for the stress- 

energy tensor for general curvature coupling in agreement with [10]. We analyzed the null 
energy condition for the specific case of conformal coupling and showed, analytically for 
large distances, and numerically everywhere, that it is obeyed. 

This implies that the averaged null energy condition is obeyed for conformal coupling in 
this system. Since the ANEC is independent of the curvature coupling, it is obeyed for any 
coupling, in accord with our conjecture in |j| that the ANEC is satisfied for all geodesies 
which pass outside any localized background potential. 

We have also given a derivation of the standard stress-energy components for a flat 
Dirichlet plate as the close-in limit of the spherical geometry considered here. The well- 
known divergence as one approaches the surface arises in this case from a summation over 
many angular momentum modes. 
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